Parameters and pitfalls in dark energy models with time varying 

equation of state 



Rahul Biswas 1 , Benjamin D. Wandelt 1 ' 2,3 
1 Department of Physics, University of Illinois at Urbana- Champaign, 1110 W. Green Street, Urbana, 

IL 61801, USA 

2 Department of Astronomy, University of Illinois at Urbana- Champaign, 1002 W.Green Street, Urbana, 

IL 61801, USA 
3 California Institute of Technology, MC 130-33, Pasadena, 
CA 91125, USA 



ABSTRACT 



Are geometrical summaries of the CMB and LSS sufficient for estimating cosmological pa- 
rameters? And how does our choice of a dark energy model impact the current constraints on 
standard cosmological parameters? We address these questions in the context of the widely used 
CPL parametrization of a time varying equation of state w in a cosmology allowing spatial cur- 
vature. We study examples of different behavior allowed in a CPL parametrization in a phase 
diagram, and relate these to effects on the observables. We examine parameter constraints in such 
a cosmology by combining WMAP5, SDSS, SNe, HST data sets by comparing the power spectra. 
We carefully quantify the differences of these constraints to those obtained by using geometrical 
summaries for the same data sets. We find that (a) using summary parameters instead of the full 
data sets give parameter constraints that are similar, but with discernible differences, (b) due 
to degeneracies, the constraints on the standard parameters broaden significantly for the same 
data sets. In particular, we find that in the context of CPL dark energy, (i) a Harrison- Zeldovich 
spectrum cannot be ruled out at 2a levels with our current data sets, and (ii) the SNe IA, HST, 
and WMAP 5 data are not sufficient to constrain spatial curvature; we additionally require the 
SDSS DR4 data to achieve this. 



Introduction 



A 


number of observations (IPerlmutter et al. 


1999; Riess et al. 19981 


Garnavich et al. 1998: 


Knoi 


3 et al. 2003; 


Tonrv et al. 2003; Riess et al. 


2004 


: Astier et al.l 


20061: Wood-Vasev et all 2007: 


Hicken et al.l l2009f)have established that the ex- 



pansion of the universe is accelerating. The cause 
is usually attributed to a currently dominant com- 
ponent called dark energy. Current data is con- 
sistent with a standard cosmological model called 
ACDM, with dark energy in the form of a cosmo- 
logical constant A. However, dark energy might, 
in fact, be a dynamical component. Indeed, there 
exist numerous models of cosmology which pro- 
duce the observed acceleration, either by postu- 
lating the the existence of one or more otherwise 
unobserved fields, or as the effects of a departure 



of gravity from General Relativity at large scales 
that cannot be ruled out by current data. There- 
fore an important objective of current and future 
observational efforts is to study the acceleration 
of the universe in different ways and detect de- 
partures in the behavior from that expected in a 
standard ACDM model. To this end it is usual 
to parametrize dark energy as a 'fluid' with its 
equation of state (EoS), and the speed of sound in 
the fluid specified independently. Generally such 
a description encompasses a large variety of phys- 
ical models if the equation of state is assumed to 
depend on the density, and the speed of sound 
depends on both the background density and the 
wavelength of perturbation. The idea is that con- 
straining these phcnomenological parameters of a 
fluid will narrow down the class of physical models 
causing the acceleration, and in particular, study 



the differences with ACDM. 

A specific time dependent parametrization of 
the EoS, and a constant speed of sound describes 
a subclass of these phenomenological models. A 
simple example is t he CPL parametrization of 
the equation of state ( Chevallier fc PolarskilbOOlt 
Linderll2003l) 



w(a) — w + Wi(l — a) 



(1) 



of a non-interacting dark energy, which has 
been adopted by th e Dark Energy Task force 
( Albrecht et al. 20061 ) in determining the relative 
importance of future experiments studying dark 
energy. Conveniently, it includes the case of a 
constant EoS (wCDM) with (w = w,w\ = 0), 
and the ACDM model (w = -1, Wi = 0). 

Parameter constraints on dark energy parame- 
ters with different time varying equation of state, 
including the CPL parametrization have been in- 
vesti gated in the context of similar recent data 
sets Jwrightll2007tlWand[200i : lLazkoz et al.ll2008t 
Davis et al.l 120071 ) . These analysis use certain 
summary parameters (the shift parameters R and 
I for the CMB, and D parameter for the Baryon 
Acoustic Oscillations, see Sec. 13.11 and references 
therein), which are intended to capture most of 
the information in data sets. In contrast, we use 
a likelihood analysis that compares the theoretical 
power spectra to spectra inferred from data, as is 
the standard practice for ACDM models. 

In this paper, we investigate the constraints 
on the parameters of a non-flat cosmology with 
a CPL dark energy from current data sets. 
These include the W MAP five year data set 
([Dunklev et al. 2008) "or the anisotropies in the 
Cosmic Microwave Background (CMB), SDSS 
DR4 data set for Luminous Red Galaxies (LRG) 
(|Tegmark et al.ll2004l 120061). the Uni on data set 
for supernovae (jKowalski et al. 20081 ). We point 
out the differences in dark energy parameter con- 
straints computed from the two approaches. Evi- 
dently, dark energy with a time varying equation 
of state behaves very differently from ACDM cos- 
mology, where the fractional density of dark en- 
ergy is tied to acceleration at a particular redshift. 
Even in the simple and comparatively benign CPL 
model this gives rise to counter-intuitive effects 
on general parameter constraints, which we study 
here. While we discuss the behavior of generic 
time varying EoS, we restrict our calculations to 



the specific case of the CPL parametrization. This 
is an appropriate example, not only because of 
the cu rrent consensus of us ing the CPL as a stan- 
dard ( Albrecht et al. I l2006h . but also because the 
it is a fairly benign evolution. Model independent 
constraints one might hope to study will include 
much stiffer variation of the EoS, where the dy- 
namical effects we describe (see Sec. [2]) could be 
more pronounced possibly leading to larger differ- 
ences with parameter constraints computed from 
summary parameters. 

In Sec. [21 we discuss the behavior of differ- 
ent relevant regions of the CPL parameter space, 
their observable signatures and their possible con- 
sequences for parameter estimation. In Sec. [3] we 
describe the details of our method to investigate 
parameter constraints derived using summary pa- 
rameters and power spectra. The results of our 
investigations are presented in Sec. |U We discuss 
our conclusions and possible implications in Sec. [5] 

2. Characteristics of Dark Energy with 
Time Varying Equation of State 

The density of a dark energy component with a 
variable equation of state w(a) evolves with scale 
factor as pue ~ a~ 3 ( 1+u,ctf ( a '), while the pressure 
P(a) ~ (1 + w(a))pY) E (a), where 



Weff(a) 



1 



J® da(l + w(a))/a 
ln(a) 

(w + wi) + (1 - a)w 1 /ln(a) 



with the last expression being valid for a CPL EoS. 
Dark energy of current density (in units of critical 
density today) SIde contributes an amount a = 
aHltt DE {l + ?>w{a))a-^ w ^^ + V towards the ac- 
celeration of the universe at a redshift z = l/a—1, 
while its density in units of matter density grows 
as ~ a 3w ° a . Thus, the presence of dark energy at a 
particular redshift modifies the background expan- 
sion. This (along with the presence of curvature) 
affects the CMB and the matter power spectrum 
(a) geometrically by altering the angular position 
of the peaks and (b) dynamically by altering the 
magnitude of the spectrum. Dark energy is be- 
lieved to be smooth (ie. its density perturbations 
do not grow at scales smaller than the Hubble 
scale), leading to less clustering of density pertur- 
bations at a particular redshift if a significant frac- 
tion of the background density is made up of dark 



2 



energy. The gravitational potentials also evolve 
differently because the background evolves differ- 
ently from a matter dominated universe. This 
leads to weaker sourcing of the growth of pertur- 
bations, thereby suppressing the number of galax- 
ies formed at a particular redshift reducing the 
matt er power spectrum i nferred from galaxy sur- 
veys ( Doran et al. 2001bl ). The change in gravita- 
tional potentials also affects the CMB power spec- 
trum through the Integrated Sachs- Wolfe (ISW) 
effect. On a plot of the matter power spectrum or 
CMB power spectrum the geometrical effects show 
up as horizontal differences (change in scales), 
while the dynamical effects show up on the ver- 
tical differences (magnitudes). 

For a dark energy characterized by a constant 
EoS, Wes(a) — w(a) — w. To explain the acceler- 
ation observed today from SN IA data, w < — 1/3 
and is close to —1 if all data sets are considered. 
This also implies that the dark energy density (in 



units of matter density) grows as 



so that 



it is negligible at earlier times. Hence, its dy- 
namical effect on the matter power spectrum is 
small, and its effect on the CMB power spectrum 
is limited to a late time ISW effect, which only 
affects the low multipoles of the spectrum. Since 
the cosmic variance at the low multipoles is high, 
the observable imprints of a redshift independent 
equation of state are limited to the geometrical ef- 
fects. For a redshift dependent equation of state, 
w(a) and therefore w c f[(a) can be quite different 
from w(0) for smaller values of a. Hence, it is 
possible to simultaneously have the observed ac- 
celeration due to very negative values of w(a = 1), 
while the ratio of densities of dark energy to mat- 
ter ~ a 3 ™t!ff(a) rema i ns significant at early times, if 
w e g(a) w 0. Thus dark energy with redshift depen- 
dent EoS can cause acceleration today, and also be 
non-negligible at early times, thereby suppressing 
the growth of structures, and affecting the CMB 
by an Early ISW Effect. It is possible for models 
with a time dependent EoS to have regions in pa- 
rameter space which are relatively indistinguish- 
able in terms of the geometric effects, but distin- 
guishable in terms of dynamical effects. This kind 
of geometrical degeneracy is likely to be more pro- 
nounced for parametrizations which allow a strong 
variation of u>(a), as that would facilitate a quicker 
transition from an equation of state around —1 
to 0. Here, we will work out the consequences 
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Fig. 1. — The upper panel shows the evolu- 
tion of w (solid) and w c g (dashed) for the set 
(w ,wi) = (-1.65116,-3.41463) (Red), (-0.541648, 
-3.07916) (Black) (0.613754, 0.333922) (Green), (- 
1.52168, 1.5576) (Blue). The lower panel shows 
the Phase Diagram of the CPL model. The re- 
gions above the horizontal solid black line are pa- 
rameters that will eventually go over to the Phan- 
tom phase causing a Big Rip, while the regions 
to the left of the solid black vertical line are in 
the phantom phase now. The region above the 
red diagonal line have significant dark energy at 
early times. The shaded regions 'cross the phan- 
tom divide', while the regions in the South East- 
ern quadrant lead to a short burst of acceleration. 
The dark energy parameters plotted in the upper 
panel are also marked. 
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in the specific example of the widely used CPL 
parametrization in Eq. [JJ which is a fairly gentle 
variation. 

The CPL parametrization is an ad-hoc parametriza- 
tion which allows the dark energy EoS w(a) (and 
w c g(a)) to asymptote between two constant values 
Wq+Wi and Wq at early and late times. In order to 
understand the behavior of this parametrization 
and the physical models it may represent, we can 
study a phase diagram on the Wq, W\ plane. In the 
upper panel of Fig. [I] we show the evolution of the 
CPL EoS w (solid lines) and w c g (dashed). It is 
useful to study this parametrization as it captures 
general features of dynamical dark energies, and 
CPL models that correspond to (wo,wx) differ- 
ent from (—1,0) are definitely dynamical models 
distinct from a cosmological constant. From the 
asymptotic behavior of the upper panel of Fig. [TJ 
we can study different phases of the CPL EoS, 
where the evolution of dark energy can be asymp- 
totically similar to different models of dark energy. 
This is studied in the lower panel of Fig. [TJ Firstly 
we can separate the regions where the dark energy 
is in the Phantom phase with equation of state 
(w(a) < — 1) which occurs for scalar field theories 
with tachyonic instabilities, or with non canonical 
kinetic en ergy terms inspired by higher deri vative 
theories (ICaldwellll20C)2 ICarroll et alJl2003h . The 
two black lines (solid) bound the regions where 
the dark energy behaves like a phantom model 
(w(a) < —1): the parameter region left of the ver- 
tical black line is currently in the phantom phase, 
while the parameter region below the diagonal 
black line was in a phantom phase at early times. 
The region above the horizontal dot-dashed blue 
line will eventually become a phantom model, 
resulting in a 'Big Rip'. The region below the 
horizontal dot-dashed line will eventually become 
'normal' with w > — 1. The region between the 
dot-dashed blue line and the diagonal black line 
(labeled quintessence) , to the right of the vertical 
black line is the only region where w(a) > — 1 at 
all times. Thus this is the only phase which may 
be similar to models of a single, non-interacting 
stable canonical scalar field. The North West- 
ern and South Eastern quadrants (shaded yellow) 
marked out by the solid black lines exhibit the 
crossing behavior from the 'phantom phase' to 
the 'normal phase' where w crosses —1; this leads 
to instabilities in the evolution of perturbations 
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Fig. 2. — The upper panel shows the ratio of the 
densities of dark energy to the density of dark mat- 
ter in models selected from different parts of the 
phase diagram in Fig. QJ The middle panel shows 
that even though dark energy dominates at large 
redshifts for the (Big Rip, Phantom) kind of dark 
energy models (blue line), it does not accelerate 
the universe at early times. The lower panel shows 
the ratio of the total density of components (mat- 
ter, dark energy, relativistic components) in units 
of the total density assuming (wq, tui)=(— 1, 0) of 
the models selected in Fig. [TJ 
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at the crossing (see Sec. 13.21 and references therein 
for details). The South Eastern quadrant shows 
a rapid change of w(a) at recent times, but de- 
cays rapidly with redshift. The dashed red line 
marks the boundary at which w(a) — at early 
times; the regions above this line have significant 
dark energy density at early times like recom- 
bination. In order to study the distinct effects, 
we mark five points on this phase diagram. The 
brown point represents a cosmological constant 
with fl m = 0.3, fi D E = 0.7, i/ = 72km/s/Mpc. 
The other points were chosen from our chains rep- 
resenting the posterior of current CMB , HST and 
SNe data using summary parameters (see Sec. 13.11 
for details). The red filled circle is a 'Phantom', 
the black diamond is 'Burst DE', the green trian- 
gle is similar to a 'Quintessence' till now, while 
the blue star will have a 'Big Rip' and is chosen to 
have significant early dark energy. We will stick 
with this color code in discussing effects, and use 
these names to label the plots when possible. 

The effect of dark energy parameters of these 
distinctive types can be seen in Fig. [2j where we 
show the impact on the background for the pa- 
rameter points marked out in Fig. [TJ The upper 
panel shows the ratio of dark energy density to 
dark matter density as a function of redshift. We 
notice that the dark energy density decays with 
redshift, except for the model 'Big Rip', which 
had significant early dark energy (cf lower panel 
Fig. [1]). The middle panel shows the acceleration 
of the universe for dark energy with the same pa- 
rameters. Clearly, the 'Phantom' and 'Big Rip' 
models with very low values of wq exhibit a super- 
accelerated phase. It is interesting to note that 
even for the 'Big Rip' where dark energy density 
is significant at early times, the universe is decel- 
erating before a redshift of 2 because its EoS (Blue 
solid curve in the upper panel of Fig. [1]) keeps in- 
creasing. Hence for these models with significant 
early dark energy, the universe can be 'dark energy 
dominated' without accelerating, with the effect of 
dark energy on background expansion being sim- 
ilar to matter. The lower panel shows the evolu- 
tion of the total density of matter, dark energy 
and relativistic components, ptot = Pm + Pde + Pr 
with redshift, in units of the corresponding quan- 
tity assuming that dark energy was a cosmological 
constant. 

The effect of these parameters on observables is 
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Fig. 3. — Effect on observables for models with 
different values of dark energy parameters. The 
parameters are the points marked in the phase 
diagram of Fig. [U chosen from different regions 
marked on the phase diagram.: The upper panel 
shows the CMB power spectrum for the param- 
eters shown. The middle panel shows the mat- 
ter power spectrum as a function of wave-number 
in units of h~ 1 Mpc. The lower panel shows the 
change Am in the apparent magnitude with red- 
shift z. However, we use a likelihood which is an- 
alytically marginalized over the absolute magni- 
tude of the supernovae; hence we only show Am 
up to a constant for each cosmology. The horizon- 
tal dashed black lines are at ±0.15 which is the 
usual intrinsic dispersion for supernovae. 
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studied in Fig. [3] for the models chosen in Fig. [T] 
Since we are interested in comparing the effect on 
the magnitudes of the spectra due to dark energy, 
we set the amplitude of the primordial power spec- 
trum to be the same for each model. This does not 
change the value of the posterior (using Likelihood 
B, see Sec. l3.ip probability of the model plotted, as 
the posterior is independent of the amplitude. The 
upper panel shows the CMB power spectrum, the 
middle panel shows the matter power spectrum, 
and the lower panel shows the apparent magni- 
tude of the SNe IA, which is a purely geometrical 
effect. The difference in positions of the angular 
positions of peaks in both the CMB and the os- 
cillations in matter power spectrum, as also the 
differences in the redshift-magnitude plot is due 
to geometrical effects. The difference in the mag- 
nitudes of the power spectrum is due to dynamical 
effects that are not captured in summary parame- 
ters (see Sec. l3.1|) . It is interesting to note that the 
models for which the CMB power spectrum is en- 
hanced have suppressed matter power spectrum. 

The comparison of the difference of the CMB 
and Matter Power Spectra inferred from the data 
with their theoretically computed counterparts 
with the scale of expected deviations leads to con- 
straints on cosmological parameters. It is pos- 
sible that a subset of this information pertain- 
ing to specific features in the spectra is almost 
as useful in constraining parameters. For exam- 
ple, many features of CMB power spectra depend 
on the cosmological parameters through very spe- 



cific f unctions of the these parameters (|Hu et al 
2001: iHu fc Dodelsonl 120021 : iDoran et alJ l2001eJ : 



Doran fe Lillevl 2002 ). This motivates the idea of 
using CMB shift parameters to summarize the in- 
formation co ntent of the CMB anisotropies, in- 

( 19971) in the context 



troduced in Bond et al 



of forecasting for sta ndard cosmologies. Further 
work on this subject ( ElgarOv fe Multamakill2007 : 



■leai 

Wang fe Mukherjed 120071 ) suggests that one can 
summarize the information the in the spectra effi- 
ciently in terms of the summary parameters R, l a , 
which relate to the position of the first peak of 
the CMB spectrum, and the spacing of the peaks 
due to acoustic oscillations. Recently such sum- 
mary parameters have become popular in study- 
ing the constraints on dark energy parameters. If 
the effect of dark energy is mostly geometrical, it 
is tempting to speed up the calculation by summa- 



rizing the CMB/LSS data in terms of a few geo- 
metrical summary parameters that describe these 
effects, instead of going through the time consum- 
ing process of theoretically computing the angular 
power spectrum. 

This approach of using summary par a meter s 
has been studied in Wang fc M ukheriee (2007); 



Elgar0v fc Multamakil(|2007l ): lDoran fc Lillevl (l2002f) 
and th e procedure nicely outlined in iKomatsu et al 
(|2008l ). One finds samples of the posterior distri- 
bution of the data set concerned for cosmological 
models with a specific form (ACDM) of dark en- 
ergy. One then uses these samples to estimate 
the set of summary parameters over the poste- 
rior distribution as well as the covariance over 
these summary parameters. IWang fc Mukheriee 



(|2007l ) also show that the summary parameters 
are comparatively weakly correlated to the other 
cosmological parameters. One then approximates 
the likelihood of the CMB data set, by a Gaussian 
distribution over these summary parameters It can 
be seen that there are three crucial assumptions 
in this procedure: 

1. The Likelihood of the CMB spectra are well 
described by CMB summary parameters. 

2. The mean and covariance of the summary 
parameters for a particular data set do not 
change appreciably when the model space is 
enlarged, and do not develop correlations. 

3. The summary parameters are relatively 
weakly correlated with the other cosmologi- 
cal parameters. 

We examine these by studying the differences in 
the parameter constraints in comparison of the 
power spectra and the use of the summary pa- 
rameters. 

3. Method 

We study parameter constraints using a suit- 
ably modified version of the Mark ov Chain Monte 
Carlo (MCMC) engine CosmoMC (|Lewis fc Bridle! 
|2002| ) to explore the Bayesian posteriors for the 
cosmological parameters. We assume an isotropic 
and homogeneous universe with dynamics dic- 
tated by standard general relativity with the den- 
sities of the background components to be deter- 
mined, take the non-baryonic dark matter to be 
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Table 1: Parameters used in the MCMC: Cosmological parameters used and the lower and upper limit on 
the flat priors on these parameters. For a wCDM model, w\ is fixed to be 0. 







e 


T 


n k 


w 


Wi 


n s 


log(10 w A s ) 


(0.005,1) 


(0.01,0.99) 


(0.5,10) 


(0.01,0.6) 


(-0.2,0.2) 


(-3,1.5) 


(-7.0,4.5) 


(0.5,1.5) 


(2.7,4) 



entirely cold and neglect effects of neutrino mass, 
assuming AT eff = 3.04 species of massless neutri- 
nos. The current background density of photons 
is also assumed to be fixed, by neglecting any 
error-bars associated with the measurement of 
the CMB temperature. We then explore possible 
values for the background densities of other com- 
ponents (baryons, cold dark matter, curvature) 
with broad, flat priors over Wb,w c ,J\, where fij 
is the density of the component in units of the 
current critical density, and u>i = h 2 fl- l7 with the 
Hubble constant H = 100/i Km/s/Mpc. The pri- 
mordial perturbations are assumed to be adiabatic 
and Gaussian distributed with a power spectrum 
Pj(fc) = A j (fc//cp)("j( fc )- 1 ). We neglect the effect of 
tensor perturbations setting the ratio r of its am- 
plitude A t to that of the scalar perturbations A s 
to be zero, and nt = 0. The scalar spectral index 
n s is assumed to be scale independent (running 
of spectral index) n T — and we choose a pivot 
scale k p = 0.05/Mpc. We also ignore the effect 
of the the SZ amplitude. We allow for a single 
re-ionization taking place at an optical depth of 
r, and parametrizing the angle subtended by 
the sound horizon at the surface of last scatter- 
ing instead of the Hubble constant Hq. There- 
fore, our model space has the fixed values u> v = 
n t — n r — r = 0, and priors on the parameters 
{ft b h 2 ,n c h 2 ,6, t, fl k , w ,wi : n s , LogflO^Ag)} sum- 
marized in Table [T] The Markov chains are as- 
sumed to have converged when the R — 1 statistic 
had been below 0.03 for a few tens of thousand 
chain steps; this results in a final R — 1 statistic 
of about - 4 x 10~ 3 - 1 x 10~ 2 . 

3.1. Data and Likelihoods 

We summarize our usage of different data sets 
and likelihoods in Tableland explain it in detail 
below. 

CMB data: We use the WMAP 5 year data in two 
different ways: 



(A) By using the publicly a vailable WMAP like- 
lihood code ( v ersion 3) (iDunklev et"aL 2008; 
iHinshaw et al. l2008t iNolta et all l2008f ). which 
compares the observation to a our theoretical com- 
putation of the CMB power spectrum, and 

(B) By using a Gaussian likelihood in the sum- 
m ary parameters Ua , R , z*l as recommended 
bv lKomatsu et alj (|2008l ) 



l f 



(1 + «*) 



-kD a (z+) 



R = sjn m H%{l + z*)D A { Zit ) 



(2) 



where z± is the rcdshift of the surface of the 
last scattering comput ed from the fitting formula 
(|Hu fc Sugivamal Il996h in terms of only densities 
of baryons (cob) and matter (tu m ) with the mean of 
the distribution taken to be the maximum Like- 
lihood values in Table 10, and the inverse covari- 
ance matrix in Table 1 1 of iKomatsu et alJ (J2008J) . 
There are slight differences in the literature about 
how the parameters l a ,R are best defined, and 
we adopt the definition of IKomatsu et al. (2008) 
as we use their numerical values for the likelihood. 

Galaxy Power Spectrum Data: We use the galaxy 
power spectrum data from the Luminous Red 
Galaxy (LRG) sample of the Sloan Digital Sky 
Survey DR4 (SDSS) 

(A) By using a modified version of the pub- 
licly available likelihood code which compares the 
matter power spectrum inferred from the data 
with the theoretical computation after analytic 
marg inalization over a sc ale independent linear 
bias ( Tegmark et al. 20061 ). We modified the code 
in order to recompute the geometric scaling at the 
redshift of the sample related to 



D v = {{l + z) 2 D\{z)cz/H{z)) 



1/3 



for the CPL model, and also to compute the 
growth function to account for the scale inde- 
pendent change in the matter power spectrum at 
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Tabic 2: Definitions of data sets I and II and Likelihoods A and B 



Data Set I 


Data Set II 


Likelihood A 


Likelihood B 


WMAP5 + SNE + HST 


WMAP5 + SNE + HST 
+ SDSS + BBN 


Summary Parameters 
Ia, Ra, z*, Da 


Power Spectra 



the mean redshift 0.35 of the LRG sample from 
the matter power spectrum at a redshift of z = 0. 
(B) By using only the geometric distance mea- 
sure r s (zd)/Dv{z), where r s is the sound hori- 
zon at the redshift of drag epoch Zd where the 
baryons were re l eased from the photons. Following 
Komatsu et al. (2008), the redshift za is computed 
through a fitting function ( Hu fc Eisensteirj|l998 ) 
which again depends only on and u> m . As rec- 
ommended, we use a Gaussian Likelihood with a 
mean of 0.1094 and a standard deviation of 0.0033. 



Supe 



iae Data: W e use the Union data set 
()Kowalski et al.l [2008) which is a compilation 
of 307 supernovae IA discovered in different 
surveys. This combines high redshift super- 
novae by the ESSENCE, SNLS and HST Goods 
surv eys, with low red s hift ones (z rj 0.02 



0.1 s ) lAstier et al.l (l2006MPerlmutter et all (|l999h : 



Riess et al.l ( 20041 ) using a weighting scheme to 



take into account the heterogeneous sources. The 
Likelihood is Gaussian distributed in the magni- 
tude space and includes covariance contributions 
due to systematic effects. However, we ignore 
lensing of supernovae. 

HST: We also incorporate the results of the Hub- 
ble Space Telescope Survey measurements of the 



Hubble constant (jFreedman et al.l 120011 ) as a 



Gaussian prior on the value of the Hubble con- 
stant H = 72± 8km/s/Mpc. 
Big Bang Nucleosynthesis: The primordial abun- 
dance of light nuclei, determined at the time of 
Big Bang Nucleosynthesis (BBN) depends on 
the baryon to photon ratio, as well as the ex- 
pansion rate at th is time ( Amsler et al. 2008t 



Cvburt et al.1 120021 ). Since, in time varying equa- 
tion of state dark energy models, dark energy 
can be non-negligible around the time of BBN 
(z ~ 10 7 ), we need to examine the effect of BBN 
on both these parameters. Among the abundances 
of light elements, the abundance of primordial He- 
lium is most sensitive to the expansion rate dur- 



ing BBN; the abundance of primordial Deuterium, 
while almost insensitive to the expansion rate, is 
extremely sensitive to the baryon to photon ra- 
tio. We ignore the abundance of 7 Li even though 
it is extremely sensitive to the baryon photon 
ratio, because of the controversies regarding sys- 
tematic uncertainties in determining the observed 
abundances of 7 Li, stemming from uncertainties 
in measurement of effective temperature of stars, 
or an unaccounted correlation betw een the metal- 
licity and estimated abundance (jCvburt et al. 
2008). We constrain the dark energy density at 



the time of BBN by using the fit equations for 
the primord ial abundance of Deuteri um, and He- 
lium from (jSimha fc Steigman 2008) in terms of 
the baryon-to-photon ratio ryio , and the ratio S of 
the Hubble parameter at times of BBN (z ~ 10 7 ) 
with a Hubble parameter for a universe completely 
dominated by relativistic degrees of freedom (pho- 
tons and massless neutrinos). The parameter rjio 
depends on the mass fractions of the light nuclei, 
but the dependence is extremely weak. We adopt 
the values of 7710 = 273. 9u>b for our Likelihood cal- 
culations. We ignore the theoretical errors in the 
fitting functions, and write a Gaussian likelihood 
using the error estimates of observed abundances 
for De and 4 He. We also note that since the 4 He 
fraction monotonically increases with time, the 
lowest values detected are a robust upper bound 
to the Helium fraction. We impose a much weaker 
hard prior on the equation of state parameters, 
as will be described and further justified in the 
Sec. H 

3.2. Theoretical Computation of Power 
Spectra 

The Likelihoods (A) for the CMB and Mat- 
ter Power Spectra data require theoretically com- 
puted values of these power spectra. We eval- 
uate these by modifying the background expan- 
sion, and the perturbation equ ations for this Dar k 
Energy model using CAMB (jLewis et all [20001. 



S 



which uses RECF AST (|Seager et al.lll999l . feoOO: 



I Wong et al. 2007) to compute the recombination 



history. The CPL parametrization allows for the 
case where the dark energy equation of state is 
less than — 1, and allows crossing of —1 dur- 
ing the evolution. For a single non-interacting 
scalar field, w > — 1. However, lower values of 
w occur for scalar fields with non-canonical ki- 



more than one field (Hu 2005: 


Huev & Wandelt 


2006 


; Zhang et al. 2006: Li et al. 


20051: Fang et al. 



2008). Thus, we follow the standard practice and 
do not rule out values of w(a) < —1. ft is known, 
that in the fluid model perturbation theory of dark 
energy, there is a runaway problem associated with 
the models which a llow crossing of —1. We follow 
the p rescription of ( Huevl[2004 : Caldwell k. Doran 
2005), which essentially assumes that this is a 
numerical artifact. While, a negative dark en- 
ergy sound speed would lead to large clustering of 
structure that is unobserved in data, we find that 
current data is insufficient to put any meaning- 
ful constraints in the open interval (0,1) (in nat- 
ural units) on the speed of sound c s for dark en- 
ergy. Therefore, we fix the value of the speed of 
sound to a reasonable choi ce of c 2 a = 1 as would 
be the case for quintesse nce ( Weller fc Lewisll2003t 
Caldwell fc Doranll2005h . 



4. Results 

We study the posterior distributions on the 
cosmological parameters and the joint posterior 
distributions of pairs of cosmological parameters. 
Use of the maximal data set available today ex- 
ploits the complementarity of different probes to 
obtain the tightest constraints on the parameters. 
On the other hand, it is useful to examine mean- 
ingful constraints from different subsets of the 
maximal data set. Constraints from subsets allow 
one to check for self-consistency of cosmological 
models since different data sets actually constrain 
different aspects of physics; and also the sepa- 
rate the constraints from assumptions inherent to 
different data sets. We use two sets of data. (I) 
WMAP 5 yr data, the Supernovae Union data set, 
the HST constraints on the Hubble constant, and 
(H) additionally BBN constraints, and the SDSS 
LRG power spectrum. The data combination (II) 
represents the maximal data set we use here. We 



present the constraints obtained for our maximal 
data set in Fig. [4] 



4.1. Features of the Posterior Distribution 

Focusing attention on the two dimensional joint 
posterior on w , W\, the left panel of Fig. [5] shows 
the posteriors of data set (I), while the right panel 
shows the posteriors for data set (II). The black 
(solid) contours are computed by the use of likeli- 
hood A, while the blue contours are computed us- 
ing Likelihood B. T he blue conto urs may be com- 
pared to Fig. 1 of IWand (j2008h . We note that 
the posterior distribution is highly non-Gaussian. 
Firstly, the posterior contours are fairly banana- 
shaped rather than ellipsoidal without the SDSS 
data. With the addition of the SDSS data, the 
elongated banana shape gets pinched off. How- 
ever, the distribution is still elongated, and even 
in the high posterior region the posterior peak 
does not seem to be well-centered. Further, the 
plots show that the posterior falls abruptly around 
the blue (dot-dashed) line w\ — ~wq. The diago- 
nal (solid) black line represents the 2a limit from 
BBN constraints; wo,w\ values lying above the 
black line are unlikely due to the BBN constraints. 
The red diagonal line is a hard prior we used to 
limit the exploration. The use of both Likelihoods 
(A) and (B) result in a dramatically lower num- 
ber of points beyond the blue (dot-dashed) line 
wo + Wi =0, though the change in B is less sharp. 
This edge in the likelihood, is entirely due to the 
CMB data and represents the edge of parameter 
space beyond which the 'dark energy' dominates 
at the redshift of the surface of last scattering, 
rather than being related to the BBN data. We 
can see this better in a histogram of asymptoti- 
cally early equation of state in Fig. [6] It shows the 
drop in the posterior sample for both Likelihoods 
B as well as A for data set I (Upper Panel) and 
data set II (Lower Panel), even though data set 
Idoes not include the BBN constraints. The lack 
of points above the (dot-dashed) blue line, further 
justifies our use of the hard prior, since the pos- 
terior is well disconnected from the hard prior. In 
our contour plot, we therefore could cut off the 
posterior contours that we would have drawn by 
smoothing the posterior densities with a Gaussian 
kernel. This sharp edge in the 2D joint posteriors 
is unlikely to go away with the addition of further 
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Fig. 4. — Marginalized one dimensional and two dimensional joint posterior distributions of our maximal 
data set (WMAP5, SDSS, HST, SNe , and BBN) on the cosmological parameters. The solid lines show the 
posteriors, while the dotted lines show the mean likelihoods. 




Fig. 5.— 2D Joint Posterior distributions given (left) data set I (WMAP5, SNe, HST) and (right) data II 
(WMAP5, SDSS, SNe, BBN, HST). The black (solid) lines were computed using the Likelihoods A comparing 
spectra, while the blue (dash-dot) lines were computed using summary parameters B. 
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Fig. 6. — Binned one dimensional posterior on 
the dark energy EoS at very early times w(0) — > 
Wo + w\ marginalized over all other parameters 
showing the sudden drop in the posterior for val- 
ues of w(0) + w(l) — 0. The upper panel shows 
the posterior given the combination data set (I), 
while the lower panel shows the posterior given the 
combination of data sets (II). The blue histogram 
shows the distribution due to the likelihood (A) 
while the magenta histogram is computed with the 
use of the likelihood (B). 



lr 

0.9 
0.8 
0.7 
0.6 

£ 0-5 
CL 

0.4 
0.3 
0.2 - 
0.1 

- 

lp 
0.9 
0.8 
0.7 
0.6 

i ° - 5 

Ok 

0.4 

0.3 

0.2 

0.1 

- 
-8 



-2 
W, 




Fig. 7. — Marginalized one dimensional poste- 
rior distributions on w\: The upper panel shows 
the posteriors for data I (WMAP5, SNe, HST) 
and lower panel shows the posteriors for data II 
(WMAP5, SDSS, SNe, BBN , HST). The black 
(solid) lines were computed using the Likelihoods 

(A) comparing spectra, while the blue (dash-dot) 
lines were computed using summary parameters 

(B) . 



Next, we look at the constraints on the dark 
energy parameter w\. In Fig. we show the 
one dimensional marginalized posterior distribu- 
tion on w\ computed from the data combinations 
I (left) and II (right) described above, according to 
the likelihoods A (black solid) and B (dot-dashed 
blue). The posteriors computed usi ng Likelihood 
B may be compared to the Fig. 2 of IWand (|2008l ) 
where they use likelihoods using as a param- 
eter; as expected these posteriors match to ex- 
tremely small differences that may be attributed 
to the slightly different data or numerical proce- 



dures. We note that these posteriors are asym- 
metric, a result of the fact that significant early 
dark energy is essentially ruled out by the CMB 
data. From our one dimensional distributions, we 
can see that while the approximate posterior (B) 
computed by using the summary have the same 
shape as the posteriors (A) computed by compar- 
ison of power spectra, there are significant differ- 
ences between likelihoods A and B in the extent 
of the tails: the posteriors using the power spec- 
tra are sharper and narrower than the posteriors 
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using the summary parameters. For our maximal 
data set (II), this translates to the tails (computed 
by using (B)) extending about twice as much as 
the tails using summary parameters (A) . From the 
one-dimensional distributions, we see that there is 
a sharp edge in the distribution of w\ at the higher 
tail, while low w\ values are are allowed. The edge 
at the tail is related to the edge in the joint pos- 
terior. The use of summary parameters using the 
likelihood (B) has a similar effect, except the dis- 
tributions are broader, and the distinction more 
significant in the low w\ tails. 

From both the one dimensional and two dimen- 
sional posteriors, we note (I) allow fairly low values 
of w\ at levels, which are ruled out by the SDSS 
data (II), or if the model is constrained to be flat. 
This is because (I) allows models with low values 
of wi 7 Hq, and flk, which are ruled out by the si- 
multaneous use of the SDSS data. This can be 
seen by studying the correlations of wq , w\ shown 
in Fig. El 

4.2. Comparison of Power Spectra and 
Summary Parameters 

The differences between the posterior distri- 
butions due to the use of different likelihoods A 
and B, are most clearly studied in a binned (un- 
smoothed) density plot over all the chains, since 
various differences can arise due to smoothing pre- 
scriptions inherent in making contours. We choose 
to study the joint posterior in the wq,w\ plane 
due to its im portance in classifyin g dark energy 
experiments ( Albrecht et alJ 2006). We present 
the differences as a density plots for both data set 
combinations I (top panel) and II (bottom panel) 
in Fig. [9j We can see that the differences are most 
appreciable for models, where the equation of state 
is close to (like matter) at asymptotically early 
times, and for the tail, where w\ has low values. 
The density plots on the left due to likelihood A 
are much sharper and narrower than the corre- 
sponding plots using B on the right. 

We also study the marginalized one dimen- 
sional posteriors in Fig. [10] on all other param- 
eters. From Fig. [7] we see that with the maxi- 
mal data set, the posterior due to Likelihood B 
(black) extends to about twice the tail of Like- 
lihood A (blue) in w\. From Fig. [101 we show 
the 1 dimensional marginalized distributions on 
all other parameters from the comparison of spec- 
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Fig. 8. — Correlations of w\ with background pa- 
rameters: 2D Joint posteriors of w\ values with Qk 
(upper panel) and H (lower panel) showing the 
68 (green and red) and 95 (yellow and blue) per- 
cent contours for data sets I (WMAP5, SNe, HST) 
and II (WMAP5, SNe, SDSS, HST, BBN) using 
the Likelihood A comparing the power spectra 

tra (ie. by using likelihood A) (blue) and the sum- 
mary parameters (red) (ie. by using likelihood B). 
We see that the posterior distributions for wo,8 
are well approximated by the summary parame- 
ters. On the other hand, comparing the blue and 
red curves in Fig. [101 we see that the posteriors 
on all other quantities (f2 m , fi^, Qde, ^b, Ho) are 
shifted significantly in comparison to the width 
of the distribution when the summary parameters 
used. As expected, the parameters related to the 
optical depth of re-ionization, scalar spectral in- 
dex, and the amplitudes of the scalar fluctuations 
(t, n s ,A s ) are not constrained at all by the use 



12 



Fig. 9. — Density Plots showing the Joint Posterior Distribution on wq,wi given the data set I 
panels), and data set II (lower panels) for likelihoods A (left) and B (right). 
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of summary parameters resulting in a distribution 
about as broad as the prior, while the spectra con- 
strain them reasonably well. This is not surprising 
as all of them affect the shape and normalization 
of the CMB power spectrum without affecting the 
distance to the surface of last scattering, or the 
angular position of the peaks. 

4.3. Impact of Relaxing the Dark Energy 
Parametrization 

In general, when a model is relaxed to a more 
general model, one expects that the constraints on 
the model parameters could become broader. This 
happens, when the model parameters are corre- 
lated to the parameters that are allowed to float in 
the more general model, but were fixed in the spe- 
cial model. In our example, ACDM and wCDM 
are special cases of the CPL model. Hence, we 
look at how the constraints get broader when w\ 
is allowed to vary. From Fig. [TQl we note the dif- 
ferences in posteriors between a wCDM model 
(black) and a CPL model (blue). We see that 
for parameters like (6, u> c ,t, A s ) and derived pa- 
rameters like (f2 m , Hq), the constraints from both 
a wCDM model and a CPL model are similar. 
However, for the parameters (wq, Q.k, ^>b, n s ) and 
derived parameters (O^crs), there are differences 
between the the posteriors of the wCDM model, 
and the CPL model. Of these, only the constraints 
on wo are well approximated by the summary pa- 
rameters. We note one cannot study the impact of 
dark energy on the standard cosmological parame- 
ters using summary parameters. These differences 
are reflected in the joint two dimensional poste- 
riors of the wCDM model (blue) and the CPL 
model (black) in Fig. [Til which show that the joint 
constraints on the CPL model are significantly dif- 
ferent. 

Of these parameters, it is interesting to note 
that the values of Q,k and n s allowed by the CPL 
model. Both these parameters are important the- 
oretically as they are strongly related to parame- 
ters in inflation. The tail of n s using the data set 
(II) are usually considered outside of the 2er limits. 
Similarly the values of shown in Fig. [5] allowed 
by the data set (I) are usually considered ruled out 
by the CMB alone. The usual statements are on 
the basis of ACDM or wCDM models, and this 
shows that relaxing constraints on the dark energy 
models can over-ride some of our intuition. The 



impact of the SDSS and BBN constraints again 
underline the importance of complementary data 
sets in this context. 

5. Summary and Discussion 

Dark energy alters the background expansion 
of the universe leading to geometric effects on the 
CMB and matter power spectrum. It also changes 
the growth of gravitational potentials leading to 
dynamical effects that modify the magnitudes of 
the power spectrum. For a constant EoS dark 
energy, the observed acceleration requires rapid 
decay of dark energy density with redshift, pre- 
cluding dynamical effects. However while dynam- 
ical effects arc relatively unimportant for constant 
EoS dark energy, they may be important for dark 
energy with time varying EoS. Hence, for such 
models, the comparison of power spectra may be 
able to distinguish between regions of parameter 
space that are degenerate to dynamical effects- 
blind summary parameters used in previous stud- 
ies. 

We study this by estimating parameter con- 
straints in a CPL cosmology using two combina- 
tions of data sets I and II (see Tabled]) WMAP five 
year data, SDSS LRG data, the Union compilation 
of supernovae data, the results of HST, and using 
the fits to BBN constraints. We do this in two dif- 
ferent methods (see Table [2J (A) where we com- 
pare the observed CMB and matter power spectra 
to our theoretical computations of these quantities 
in a cosmology with a CPL parametrized dark en- 
ergy using modified versions of the publicly avail- 
able WMAP and SDSS Likelihood codes, and (B) 
where we use a Gaussian likelihood in the three 
summary parameters for CMB, and the BAO sum- 
mary parameters as reported in previous studies 
in the literature. 

Differences between Constraints : Qualitative fea- 
tures in the wq,wi joint posteriors are similar 
when either likelihood (A) or likelihood (B) are 
used, though there are quantitative differences. 
The differences in the wq, w\ joint posterior which 
is used to compute the dark energy Figure of Merit 
recommended by the Dark Energy Task force to 
rank the importance of experiments are studied 
in Fig. [5J This is an un-smoothed density plot to 
clearly show the differences between using Likeli- 
hoods A and B. Clearly, likelihood B using sum- 
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Fig. 10. — Marginalized one dimensional constraints on cosmological parameters given all current data sets 
(II) (WMAP5, SDSS, SNe, HST, BBN) using Likelihoods (A) for a wCDM model (solid, black), (A) for a 
CPL model (dashed, blue), and (B) for a CPL model (dotted, red) 
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mary parameters is good as an approximate likeli- 
hood for the CPL model for the purpose of study- 
ing dark energy parameters; however the likeli- 
hood (A) gives sharper and narrower posterior dis- 
tributions. We study the differences in one dimen- 
sional marginalized constraints on the cosmologi- 
cal parameters when one uses the Likelihoods A 
and B in Fig. [10] We find that the distributions 
due to B are mostly broader, but also include sig- 
nificant biases in some cases, further Likelihood 
B does not constrain the parameters related to 
n Sl T,A s at all. It is not surprising that the likeli- 
hood comparing power spectra is more informative 
than one using summary parameters, as the latter 
only use a subset of the information available in 
the former. 

Features of the Posterior Distributions of the CPL 
model : Using the method (A), we show the fea- 
tures of the current constraints on the CPL dark 
energy parameters wq, and wi in Fig. O The 
main features are (i) a sharp drop in the poste- 
rior for models that allow significant early dark 
energy demonstrated in Fig. [6] and (ii) a long tail 
for the combination of data sets (I) in the direc- 
tion of low wi values, for which the dark energy 
equation of state changes decreases in the past, 
resulting in a specific fractional density of dark 
energy causing more acceleration in the past; but 
also implying that dark energy itself decays away 
more rapidly, and (iii) the distribution is fairly 
non-Gaussian. The drop in the posterior for pa- 
rameters, described in (i) that allow for signifi- 
cant dark energy happens for a combination of 
geometrical and dynamical effects. Dynamical ef- 
fects come from an Early ISW effect, and a dif- 
ferent sourcing of the growth of matter pertur- 
bations due to different growth behavior of the 
potential. In Fig. [T] dark energy parameters of 
the kind described in (ii) were called "burst dark 
energy" since they lead to short burst of accelera- 
tion as shown in Fig. [2l Studying the correlations 
of the dark energy parameters with other back- 
ground parameters in Fig. suggests that such 
models are allowed if the flatness and Hubble con- 
stant values are also low, which are almost ruled 
out when we include the SDSS data. The non- 
Gaussian banana-shaped posterior on wo,wi due 
to the data sets (I) are pinched off when data set 
(II) is used, however the distribution is still not 
very ellipsoidal. This is typical of posteriors of 



ill-constrained parameters, and can change with 
the addition of higher quality data. However, in 
this case the distribution is unlikely to be ellip- 
soidal if it extends to early dark energy cutoff de- 
scribed in (i), which is not too far from the peak 
of the distribution. These features are inherited in 
the marginalized one dimensional posterior of w\ , 
where we see an asymmetric distribution with a 
sharp drop in the distribution at high values of w\ 
and a long tail in the low values of w± in Fig. [7J 
Impact on Standard Cosmological Parameters : 
We study the differences in one dimensional 
marginalized constraints on the parameters in a 
wCDM model, when the model is relaxed to a 
CPL model in Fig. [TUJ We find that the con- 
straints on Wo, n Sl (jg, cob are different in a CPL 
model (blue) from their counterparts in a wCDM 
model (Black). We also compare the two dimen- 
sional joint posteriors of these cosmological pa- 
rameters for a wCDM model and a CPL model 
in Fig. 1 1 11 where the posteriors also change when 
these parameters are involved. In particular, we 
note from Fig. [10] that in a CPL model, there is a 
tail on the higher side of the posterior distributions 
of our maximal data set. Values above unity are 
allowed, in contrast to the situation for a KCDM 
model. In Fig. we show that, in contrast to 
the situation in KCDM model with spatial cur- 
vature, where data set I (WMAP5, HST, SNe) 
constrains the flatness parameters, this is not pos- 
sible in CPL models; Addition of the SDSS and 
BBN constraints are crucial to pinching off the 
banana in Fig. [8] As discussed, the posteriors on 
these parameters using the summary are not very 
good; so these questions cannot be addressed by 
summary parameters. 

While the analysis presented is for the CPL 
paramctrization, one should remember that EoS of 
dark energy might be quite different. We note that 
the dynamical effects are likely to be more pro- 
nounced for equations of state that have a stiffer 
evolution with redshift. On the other hand, if 
the EoS has a functional form significantly differ- 
ent from a CPL paramet rization, the compute d 
constraints may be biased (jLinden fc Virev 2008). 
Attempts to circumvent this problem have been 
made in terms of increasing the number of param- 
eters describing EoS to enlarge the model space 
further, with an ultimate goal of making the de- 
scription 'model independent'. This number of 



17 



parameters that can be added is limited by the 
absence of tracers of cosmic evolutio n at high 



redsh ifts (jLinder &; Hutererl 120051 : ISarkar et al 
2008) even with the addition of futuristic SNe 



data. Typically, these analyzes are made pos- 
sible by the degeneracies between other back- 
ground parameters and the EoS by either invoking 
CMB constraints at very hig h redshift in terms 
of the summary paramet ers (jSarkar et al. 2008; 



IWang fc Mukherieel2007l) or by taking the cosmol 



ogy to be flat (IFav fc Tavakoll2006t lAlam et al 
l2004tlShafieloo et al.ll2006h . While our results sug 
gest that such a use of summary parameters may 
be safe for CPL like parametrizations, it is un- 
clear how good they are for other models intended 
to be included in these enlarged sets. While the 
assumption of flatness can be justified on the ba- 
sis of an inflation prior, we show that it cannot be 
justified using the data. The WMAP5 constraints 
on flatness for a universe with a cosmological con- 
stant significantly worsen in a CPL model and one 
intuitively expects them to be still broader for a 
dark energy with a parameter independent EoS. 

In comparing the matter power spectrum, 
we use analytic marginalization over a linear 
scale independent bias; thus parameter estima- 
tion from the matter power spectrum might be 
biased if there is really a scale dependent bias 
(jRassat et alJl2008l ) . It is clear from their analysis 
as well as our results, that there is information 
in the power spectrum that is not encoded in the 
Baryon Acoustic Oscillations. Hence, it is yet 
another reason to study the bias of galaxies, so 
that one may be able to extract the maximum 
information from galaxy surveys. 

While we have showed that a likelihood com- 
paring power spectra is more informative that like- 
lihoods comparing summary parameters even for 
CPL models, the computation of power spectra 
involves running computationally intensive Boltz- 
mann code repeatedly to obtain a large number 
of posterior samples, while the computation of 
summary parameters is extremely fast. To give 
a quantitative idea, it takes about six hours to 
get a thousand chain steps on a single proces- 
sor using Likelihood A, while it takes about six 
minutes to do the same using Likelihood B. Since 
such computations are quite doable for a particu- 
lar model, one should use a comparison of power 
spectra to do precision cosmology, while summary 



parameters can still be used to explore new mod- 
els and get rough estimates. It might be possible 
to make the likelihood using summary parameters 
more informative by adding further summary pa- 
rameters; an example of th is approach has been 
followed in (|Wrighdl2007h . where several addi- 
tional summary parameters have been used. How- 
ever, if it might be expected that a particular 
parametrization (such as the CPL model today) 
will have repeated use, one can also train B oltz- 



mann accelerators (|Fendt k Wandeldl2007al lbT) to 
efficiently and accurately compute the CMB and 
matter power spectrum. If the parametrization is 
studied enough, the work in training the Boltz- 
mann accelerators would be offset by the gain in 
time for computation of accurate constraints. 
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